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In light of the large amount of new experimental data, we revisit the determination of \ Vcb\ and mj, 
from inclusive semileptonic and radiative B decays. We study shape variables to order AqQ^/mj and 
as/3o, and include the order as AQco/mi, correction to the hadron mass spectrum in semileptonic 
decay, which improves the agreement with the data. We focus on the IS and kinetic mass schemes 
for the 6 quark, with and without expanding mt — nic in HQET. We perform fits to all available 
data from BABAR, BELLE, CDF, CLEO, and DELPHI, discuss the theoretical uncertainties, and 
compare with earlier results. We find \Vcb\ = (41.4 ± 0.6 ± Q.lrg ) x 10"^ and ml^ = 4.68±0.03 GeV, 
including our estimate of the theoretical uncertainty in the fit. 



I. INTRODUCTION 

In the last few years there has been intense theoret- 
ical and experimental activity directed toward a pre- 
cise determination of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix element \Vcb\ from combined fits to inclu- 
sive semileptonic B decay distributions [1-4]. The idea 
is that using the operator product expansion (OPE), suf- 
ficiently inclusive observables can be predicted in terms 
of iV^cfcl, the b quark mass, m?,, and a few nonperturba- 
tive matrix elements that enter at order AqQ^/m^ and 
higher orders. One then extracts these parameters and 
I Vet I from shapes of B decay spectra and the semileptonic 
B decay rate. This program also tests the consistency of 
the theory and the accuracy of the theoretical predictions 
for inclusive decay rates. This is important also for the 
determination of \Vub\i whose error is a major uncertainty 
in the overall constraints on the unitarity triangle. 

The OPE shows that in the mf, ^ Aqcd limit inclusive 
B decay rates are equal to the h quark decay rates [5, 6], 
and the corrections are suppressed by powers of «<; and 
AqcD/^fc- High-precision comparison of theory and ex- 
periment requires a precise determination of the heavy 
quark masses, as well as the nonperturbative matrix el- 
ements that enter the expansion. These are Ai_2, which 
parameterize the nonperturbative corrections to inclu- 
sive observables at ©(Aqq^/to^). At order Aqq^/to^, 
six new matrix elements occur, usually denoted by pi^2 
and Ti,2,3,4- 

In this paper, we perform a global fit to the available 
inclusive decay observables from BABAR, BELLE, CDF, 
CLEO, and DELPHI, including theoretical expressions 
computed to order a^/So, Q^s Aqcd /^Tib and Aq^p/m^. A 
potential source of uncertainty in the OPE predictions 
is the size of possible violations of quark-hadron dual- 
ity [7]. Studying B decay distributions is the best way 



to constrain these effects experimentally, since it should 
infiuence the relationship between shape variables of dif- 
ferent spectra. We find that at the current level of pre- 
cision, there is excellent agreement between theory and 
experiment, with no evidence for violations of duality in 
inclusive b c decays. 

A previous analysis of the experimental data was pre- 
sented in 2002 [1]. There has been considerable new data 
since then, which has been included in the present anal- 
ysis, and reduces the errors on |Vcb| and mi,. In addi- 
tion, the Aqcd /mt corrections to the hadronic invari- 
ant mass spectrum as a function of the lepton energy 
cut have now been computed [8], and are included in 
the present analysis. This reduces the theoretical uncer- 
tainty on the hadronic mass moments. We also compare 
our results with other recent analyses [2, 4, 9]. 



II. POSSIBLE SCHEMES 

The inclusive B decay spectra depend on the masses of 
the b and c quarks, which can be treated in many different 
ways. The b quark is treated as heavy, and theoretical 
computations for i?^*) decays are done as an expansion 
in powers of AQco/^^b- The use of the l/mt expansion 
is common to all methods. 

The decay rates for B Xc decay depend on the 
mass of the c quark, for example, through its effect on 
the decay phase space. One can treat the c quark as a 
heavy quark. This allows one to compute the D*-*-* meson 
masses as an expansion in powers of Aqcd /'tis ■ The ob- 
served masses can be used to determine m^. Since 
the computations are performed to Aq^jj/mj?, this intro- 
duces errors of fractional order Aq^q/to^ in TOc, which 
gives fractional errors of order Aqqy) / i'm-b^c) f^'^ in- 
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elusive B deeay rates, sinee eharm mass effeets first enter 
at order mfjm^. In this method, one starts with the pa- 
rameters \Vcb\, TOfc, Tic, Ai,2, Pi, 2 and Ti^^. The B* — B, 
D* — D and B — D mass differences can be used to ehmi- 
nate mb — mc, A2 and p2- Only mass differences are used 
to avoid introducing the parameter A of order Aqcd; 
thus we do not use the B meson mass to eliminate rufe. 
Three linear combinations of the four TJ's occur in in- 
clusive B decays, and the remaining linear combination 
would be needed for inclusive B* decays. In summary 
the parameters used are (i) \Vcb\; (ii) one parameter of 
order the quark mass: to;,; (iii) one parameter of order 
^QCD- '^li ^^'^ (i^) iouT parameters of order Aq^j-,: pi, 
7i — 374, 72 + 74, 73 + 374. These seven parameters are 
determined by a global fit to moments of the B decay dis- 
tributions, and the semileptonic branching ratio. This is 
the procedure used in Ref. [1]. 

An alternative approach is to avoid using the 1/mc 
expansion for the charm quark [9], since it intro- 
duces Aqcd/toc corrections, which are larger than the 
AQCD/T^6 corrections of the l/mt expansion. In this case 
heavy quark effective theory (HQET) can no longer be 
used for the c quark system, and there are no constraints 
on rric from the D and D* meson masses. At the same 
time, it is not necessary to expand heavy meson states in 
an expansion in l/mb_c, so that the time-ordered prod- 
ucts 7i_4 can be dropped. With this procedure, one has 
in addition to (i) \Vcb\; (ii) two parameters of order the 
quark mass: nib.c] (iii) two parameters of order Aqqj-,: 
Ai,2; and (iv) two parameters of order Aqqj-,: pi^2- The 
number of parameters is the same whether or not one 
expands in 1/mo. If one does not expand, two parame- 
ters of order Aq^j-, are replaced by two lower order pa- 
rameters, one of order the quark mass, and one of order 
Aqqq. The expansion parameters, such as Ai,2 are not 
the same in the two approaches. The values of Ai,2 not 
expanding the states in l/mq arc the values of Ai_2 plus 
various time-ordered products 7i — 74 when one expands 
the states in powers of l/niQ. 

In addition to the choice of expanding or not expand- 
ing in 1/mc, one also has a choice of possible quark mass 
schemes. It has long been known that a "threshold mass" 
definition for mf, is preferred over both the pole and the 
MS schemes, and it was shown in Ref. [1] that the ex- 
pansions are indeed bettor behaved in the IS" [10, 11] or 
PS [12] schemes for nib- If one expands in 1/mc, then 
rric is eliminated through use of the meson masses, and 
does not enter the final results. If one does not expand 
in 1/mc, then rric is a fit parameter. In this method, mc 
is treated as much lighter than rub, so the charm quark 
mass is chosen to be rndp-), the MS mass rcnormalized 
at a scale ~ m^. This is similar to how strange quark 
mass effects could be included in B — > Xg-y decay. In our 
computation, we will choose the scale fj, = rrib- 

In addition to the IS, PS, pole and MS schemes, we 
have also used the kinetic scheme mass for the 6-quark, 
m,|j'"(/x), rcnormalized at a low scale /i ^ 1 GcV. The 
scale n enters the definition of the kinetic mass, and 



should not be confused with the scale parameter in di- 
mensional regularization. The relation between the pole 
and kinetic masses is computed as a perturbative ex- 
pansion in powers of ««(/«), so one cannot make /x too 
small. In the kinetic scheme [9] the definitions jj^ = 
-\i + 0{as), Mg = 3A2, Pd = pi+0{as), andp^g = 3p2 
are used. 

One cannot decide which scheme is best by counting 
parameters, or by assuming that not expanding in 1/mc 
is better than expanding in 1/mc. Ultimately, what mat- 
ters is the accuracy to which experimentally measured 
quantities can be reliably computed with currently avail- 
able techniques. For example, full QCD has two parame- 
ters, mf, and rric, which can be fixed using the B and 
D meson masses. [Unfortunately, this is not possible 
to less than 1% precision at the present time.] Then 
one can predict all inclusive B decays, as well as the B* 
and D* masses with no parameters. This would be the 
"best" method to use — unfortunately, we cannot accu- 
rately compute the desired quantities reliably in QCD. 
At present, it is better to use the HQET expansion in 
1/mb and l/rric, with 6 parameters, and compute to or- 
der l/m,Q. In the (distant) future, it could well bo that 
using full QCD, with no parameters, is the best method 
to use. 

Wo have done a fit to the experimental data using 
11 schemes: the 15, PS, pole, MS and kinetic schemes 
expanding in 1/mc, not expanding in 1/mc and using 
mc(mfc), and finally, not expanding in l/rUc and using 
the kinetic scheme for both m,b and rric- In addition, 
the PS and kinetic schemes introduce a scale p, which 
is sometimes called the factorization scale. We have 
also examined the factorization scale dependence which 
is present in those two schemes. We confirm the con- 
clusions of Ref. [1], that the pole and MS schemes are 
significantly worse than the threshold mass schemes, as 
expected theoretically. This holds regardless of whether 
or not one expands in l/rric- We recommend that these 
schemes be avoided for high precision fits to inclusive B 
decays. We also find that the PS scheme gives results 
comparable to those of the IS* scheme (both expanding 
and not expanding in 1/mc), and that the PS scheme 
results do not significantly depend on the choice of fac- 
torization scale. We compared the PS scheme with the 
IS scheme in Ref. [1], and do not repeat the results here. 

Based on the above discussion, we present our results 
in five different mass schemes, using: 

1. and expand m^ — rric in terms of HQET 
matrix elements [Scheme 15exp], 

2. and rfic{mb) and no time ordered products 

[Scheme IS'no], 

3. m^™(/z = 1 GeV) and expand rrib — rric [Scheme 
kiuEXp], 

4. m|j'"(/x = 1 GcV) and mdmb) [Scheme kinNo]; 

5. m^'"(/i = 1 GeV) and m^'°(/x = 1 GeV) [Scheme 
kinuc]- (1) 
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Schemes IS'exp and kinExp contain time ordered prod- 
ucts at order Ag^-ijj/mjJ, while they are absent from 
l^No, kinNO) and kinuG- As discussed, the latter three 
schemes have the charm quark mass as an additional pa- 
rameter at leading order in Aqqd /nib. Scheme IS'exp is 
that used in Ref. [1], while scheme kinuG is that used in 
Ref. [9]. 



III. SHAPE VARIABLES AND THE DATA 

We study three different distributions, the charged lep- 
ton energy spectrum [13-16] and the hadronic invariant 
mass spectrum [8, 15, 17, 18] in semileptonic B — > Xciu 
decays, and the photon spectrum in B ^ Xg-y [19 22]. 
The theoretical predictions for these (as well as for the 
semileptonic B — > XciP rate [23]) are known to order 
a^Po and Ag^j-j/mj^, where /3q = 11 — 2ny/3 is the coef- 
ficient of the first term in the QCD /3-function. For the 
B Xii? branching rate, we use the average of the B± 
and B° data as quoted in the PDG [24], ^ 



B{B^XlD) = 10.73 ±0.28%. 



(2) 



We apply a relative —2% correction to B{B — > Xiu) to 
account for the B — > XuiP fraction, and so use 



B{B XJP) = 0.98 B{B XlV). 



(3) 



The uncertainty of \Vub\ is not a dominant error in 

B{B ^ Xctv). The fit result for \Vcb\ depends not only 
on B{B — > Xc^v), but also on the partial semileptonic 
branching ratios measured by the BABAR Collabora- 
tion [25], which have smaller errors than Eq. (2). The 
published BABAR results have already been corrected 
for B — > Xu£P contamination. 



A. Lifetime 

The value of \Vcb\ depends on the B meson lifetimes. 
The ratio of B+ and B" Ufetimes is t+/to = 1.086 ± 
0.017 [24]. Isospin violation in the B meson semileptonic 
width is expected to be smaller than both r+ /tq and the 
uncertainties in the current analysis. The w 8% isospin 
violation in the lifetimes is probably due to the nonlep- 
tonic decay channels. 

An additional source of isospin violation in the experi- 
mental measurements is through the production rates of 
B~^ and B^ mesons, which is expected to be of the order 
of a few percent to perhaps as large as 10% [26]. Let /+ 
and /o be the fraction of S+ and B^ mesons produced 



^ It would be inconsistent to use the average b hadron semileptonic 
rate (including Bs and A;, states) , since hadronic matrix elements 
have different values in the B/B* system, and in the Bg/B* or 
Ak- 



in T{4S) decay, with /+ + /o = 1. Then the measured 
semileptonic branching ratios are 



(4) 



where B and Fgi are computed with the same lepton en- 
ergy cut. In writing Eq. (4), we have used the fact that 
isospin violation in the semileptonic rates is small, so that 
the same value for Fgi is used for both and 

The measured semileptonic branching ratios can thus 
be written as 



B = Teff Fsl 

in terms of the effective lifetime 

TeS = f+T+ + /oTq. 

One can rewrite this as 

T++TO , (/+ - /o)(t+ - To) 



Teff 



+ 



(5) 



(6) 



(7) 



Using the PDG 2004 lifetime values, and the measured 
/+//o ratio [27] gives 



Teff = 1.60 ± 0.01 pS 



(8) 



where the contribution from the second term in Eq. (7) 
is negligible to both the value and the error. 



B. Lepton Moments 

For the charged lepton energy spectrum we define the 
integrals 

RniEcut, M) = / {E(- — dEt , (9) 

-'Scut 

where dT/dE^ is the spectrum in the B rest frame and 
Ecut is a lower cut on the lepton energy. Moments of the 
lepton energy spectrum with a lepton energy cut Ecut are 
given by 



Rn{Eciii, 0) 



and central moments by 



{{Ee-{Ee)r) 



RnjEcut, {Ei)E^^^) 
-Ro(-Ecut, 0) 



(10) 



(11) 



which can be determined as a linear combination of the 
non-central moments. 

The BABAR Collaboration [25] measured the partial 
branching fraction TBRo{Ecnt,^), the mean lepton en- 
ergy (Ei)^^^^, and the second and third central moments 
{{El — {Ei))")^ for n = 2,3, each for lepton energy 
cuts of iJcut = 0^6, 0.8, 1.0, 1.2, and 1.5 GeV. 
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The BELLE Collaboration [28] measured the mean lep- 
ton energy and the second central moment for lepton en- 
ergy cuts of Ecnt = 0.6, 0.8, 1.0, 1.2, and 1.5 GeV. 

The CLEO Collaboration [3, 29] measured the mean 
lepton energy and second central moment (variance) for 
E^^t = 0.7 - 1.6 GeV in steps of 0.1 GeV. 

The DELPHI Collaboration [30] measured the mean 
lepton energy, and the n = 2, 3 central moments, all with 
no energy cut. 

In total, we have 53 experimental quantities from the 
lepton moments, 20 from BABAR, 10 from BELLE, 20 
from CLEO, and 3 from DELPHI. 



C. Hadron Moments 



the expansion of (m^'^^) is wcll-behavcd. Thus, the 
calculations of the half-integer moments as presented 
in [8, 9] do not correspond to a power series in l/mb,c 
in the rric <C irib limit and omitted terms are only power 
suppressed in the small velocity limit. In addition, the 
BLM corrections to these moments arc currently un- 
known, because they require the BLM contribution from 
the virtual terms, which have not been computed. For 
these reasons we will not use these half-integer moments 
in the fit, but will compare the fit results with the mea- 
sured values. Omitting the half- integer moments, there 
are 16 data points from BABAR, 8 from BELLE, 2 from 
CDF, 4 from CLEO, and 4 from DELPHI, for a total of 
34 measurements. 



For the B 
we define 



Xclv hadronic invariant mass spectrum, 



D. Photon Spectrum 



Ie, 



dr 

dm|- 



dm 



X 



Ie, 



dr 

dm|- 



(12) 



dm. 



where £^cut is again the cut on the lepton energy. Some- 
times mj~i = [{ttid + 3to_d»)/4]^ is subtracted out in the 
definitions, {{mj^ — rfijj)"'), or the measurements of the 
normal moments are quoted, ((m^ — (m^))"), but these 
can easily be computed from (m^). 

The BABAR Collaboration [31] measured the mean 
values of mx, mj^, m\ and (i.e., n = 1/2, 1, 3/2, 2 
moments) for lepton energy cuts Scut = 0.9 — 1.6 GeV in 
steps of 0.1 GeV. 

The BELLE Collaboration [32] measured the mean 
values of mx and m\ for lepton energy cuts E'cut = 
0.9 - 1.6 GeV in steps of 0.1 GeV. 

The CDF Collaboration [33] measured the mean value 
of m\ it;s variance, with a lepton energy cut i^cut = 
0.7 GeV. 

The CLEO Collaboration [34] measured the mean 
value of m\ — rri\, and the variance of m\ lepton 
energy cuts of 1.0 and 1.5 GeV. 

The DELPHI Collaboration [30] measured the mean 
value of m\ ~ {^\: ~ ^i))^) the variance of m\, 
and the third central moment of m\, all with no energy 
cut. 

Recently half-integer moments of the m\ spectrum [8, 
9] have received some attention. While non-integer mo- 
ments of the lepton energy spectrum have been computed 
in a power series in l/mft [35], this is not true for frac- 
tional moments of the spectrum. In [8, 9] expres- 
sions for the half-integer moments were proposed which 
involve expansions that were claimed to be in powers of 
AqcD/m-c- However, in the limit TOc ^ mi, (i.e., TOc of 
order ^TO(,Aqcd or less) , the higher order terms in these 
expansion scale with powers of rrib Aqcd /'Tic' which in 
this limit is of order unity or larger. On the other hand, 
in the small velocity limit, mi, ~ nxc ^ nib—nic ^ Aqcd, 



For B —i Xgj, we define 



(-^^)Bcut = 



L 



dr 



^ dE^ 



dE^ 



dr 

d£^ 



(13) 



dE^ 



where dT/dEj is the photon spectrum in the B rest 
frame, and i^cut is the photon energy cut. In this case 
the variance, {{E^ - {E^))^) = {E^) - (E^)^, is often 
used instead of the second moment, and higher moments 
are not used as they are very sensitive to the boost of 
the B meson in the T rest frame (though this is ab- 
sent if dr/dE"^ is reconstructed from a measurement of 
dT/dEmxJ and to the details of the shape function. 
(£■") are known to order a^/So [21] and A.QCD/''^b [22]. 
These moments are expected to be described by the OPE 
once mB/2 - Ecut > Aqcd- Precisely how low i?cut 
has to be to trust the results can only be decided by 
studying the data as a function of E^ut', one may expect 
that E'cut = 1.8 GeV available at present [36] is sufficient. 
Note that the perturbative corrections included are sen- 
sitive to the mc-dependence of the b ccs four-quark 
operator (O2) contribution. This is a particularly large 
effect in the O2 — O7 interference [21], but its relative 
influence on the moments of the spectrum is less severe 
than that on the total decay rate. 

We use the BELLE [36], and CLEO [37] measurements 
of the mean photon energy and variance, with photon 
energy cuts of 1.8 and 2.0 GeV, respectively, and the 
BABAR measurement [38] of the mean photon energy 
with a cut of 2.094 GeV for a total of 5 measurements. 



IV. FIT PROCEDURE 

As discussed in Sec. II, there are many ways to treat 
the quark masses and hadronic matrix elements that oc- 
cur in the OPE results for the spectra. In the schemes 
where mi, — rric is expanded in HQET (such as 1>Sexp 
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and kinExp), the theoretical expressions for the shape 
variables defined in Eqs. (9), (12), and (13) include 17 
terms 

Xb,„, = + A + A2 + A^ + X(^) Ai 
+ A2 + X(^) AiA + X(^) AaA + X^^'^ pi 

+ T4 + e + el^M + e A , 

(14) 

while in the schemes when rric is treated as an indepen- 
dent free parameter (such as I^nOj kiuNO) and kinyc)) 
we have 22 terms 

Ye,^, = + a + f(3) a, + y(4) h? + y(5) aa^ 

+ a2 + y(7) A3 ^ y(8) a2a^ + yO) AA^ 
+ A3 + y(ii) Ai + y(i2) A2 + y(^3) AiA 

+ y(") A2A + y(i5) AiA, + y(i6) A2A, + y^^^) pi 

+ y p2 + e + y(20) el + y(2i) e A 

+ y(22)eAe. (15) 

In Eqs. (14) and (15) A and Ag are respectively the dif- 
ferences between the b and c quark masses and their ref- 
erence values about which we expand. The coefficients 
X^*^) and y^'') are functions of i^cut, and {Xe^^^, Ye^^C) 
are any of the experimental observables discussed earlier. 
The parameter e = 1 counts powers of a^. Wc have used 
Oisirub) = 0.22. The strong coupling constant is not a 
free parameter, but is determined from other measure- 
ments such as the hadronic width of the Z. The hadron 
and lepton moments are integrals of the same triple dif- 
ferential decay rate with different weighting factors. The 
use of different values of as for the hadron and lepton 
moments, as done in Ref. [9], is an ad hoc choice. 

V. THE FIT 

We use the program MINUIT to perform a global fit 
to all observables introduced in Sec. Ill in each of the 11 
schemes mentioned in Sec. II. There are a total of 92 lep- 
ton, hadron, and photon moments, plus the semileptonic 
width, to be fit using 7 parameters, so the fit has = 86 
degrees of freedom. 

To evaluate the required for the fit, wc include both 
experimental and theoretical uncertainties. For the ex- 
perimental uncertainties we use the full correlation ma- 
trix for the observables from a given differential spectrum 
as published by the experimental collaborations. In ad- 
dition to these experimental uncertainties there are the- 
oretical uncertainties, which correspond to how well we 
expect to be able to compute each observable theoreti- 
cally. For a given observable, our treatment of theoretical 
uncertainties is similar to that in Ref. [1]. 

It is important to include theoretical uncertainties in 
the fit, since not all quantities can be computed with 



the same precision. We have treated theoretical errors 
as though they have a normal distribution with zero 
mean, c uid standard deviation equal to the error esti- 
mate.^ Strictly speaking, the theoretical formula has 
some definite higher order correction, which is at present 
unknown. One can then view the normal distribution 
used for the theory value as the prior distribution in a 
Bayesian analysis. The way in which theoretical errors 
are included is a matter of choice, and there is no unique 
prescription. 

We now discuss in detail the theoretical uncertainties 
included in the fit. Those who find this procedure ab- 
horrent can skip the entire discussion, since we will also 
present results not including theory errors. 



A. Theory uncertainties 

Theoretical uncertainties in inclusive observables as 
discussed here originate from four main sources. First, 
there are uncertainties due to uncalculated power cor- 
rections. For schemes 15exp and kiuExp, these are of 
order Aq(^j-, /(m^m^) ^ 0.001, while for schemes 15no 
and kinNo where no l/rric expansion is performed, these 
are of order Aq^jj/m^ ~ 0.0001. Next, there are un- 
certainties due to uncalculated higher order perturbative 
terms. In particular, the full two loop result propor- 
tional to ag/(47r)^ ~ 0.0003 is not available. An alterna- 
tive way to estimate these perturbative uncertainties is 
by the size of the smallest term computed in the scries, 
which is the term proportional to ag/3o- We choose here 
to use half of this last computed term as an estimate of 
the uncertainty. There are also uncalculated effects of 
order {as/ '^'K)t^^QYi/m1 ~ 0.0002. Finally, there is an 
uncertainty originating from effects not included in the 
OPE in the first place. Such effects sometimes go un- 
der the name "duality violation," and are very hard to 
quantify. For this reason, we do not include an explicit 
contribution to the overall theoretical uncertainty from 
such effects. If duality violation would be larger than 
the other theoretical uncertainties they would give rise 
to a poor fit to the data. To determine the uncertain- 
tics for dimensionful quantities such as the moments con- 
sidered here, we have to multiply these numbers by the 
appropriate dimensionful quantity. This number is ob- 
tained from dimensional analysis, and we use for the n'th 
hadronic moment (m^)"/"' while we use (mB/2)"/„ for 
the n'th leptonic moment. The factors /„ are chosen to 



^ This is the same procedure as that used in doing a fit to the 
fundamental constants [39] . An example which maJces clear why 
theoretical errors should be included is: The Hydrogen hyperfine 
splitting is measured to 14 digits, but has only been computed 
to 7 digits. The Positronium hyperfine splitting is measured 
and computed to 8 digits. It would not be proper to give the 
H hyperfine splitting a weight 10® larger than the Ps hyperfine 
splitting in a global fit to the fundamental constants. 
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be /o = /i = 1, /2 = 1/4 and ,h = l/(6\/3). The values 
for /2 and /a are the maximum allowed values for the sec- 
ond and third central moments (variance and skewness) 
for a probability distribution on the interval [0, 1]. 

The complete BLM piece has not been computed for 
the non-integer hadronic moments. The perturbative un- 
certainty is therefore dominated by this contribution of 
order /3oa^/(47r)2 ~ 0.003. We will use A ~ 0.003 for the 
non-integer hadronic moments when we compare experi- 
ment with theory. 

For the hadronic mass and lepton energy moments, 
which depend on the value of the cut on the lepton en- 
ergy, we have to decide how to treat the correlation of the 
theoretical uncertainties. In the global fit by the BABAR 
Collaboration [4] , the theory errors for a given observable 
with different cuts on Ei were treated as 100% correlated. 
This ignores the fact that the higher order terms omitted 
in the OPE depend on the lepton energy cut. In Ref. [1], 
only the two extreme values of the lepton energy cut were 
included in the fit, and the correlation of the theory un- 
certainties was neglected. Here we take the correlation 
of the theoretical uncertainties to be given by the cor- 
relation between the experimental measurements, which 
captures the correlations due to the fact that observables 
with different cuts share some common events. 

For the photon energy moments an additional source 
of uncertainty is the fact that the presence of any exper- 
imentally sensible value for Ecut affects the mean pho- 
ton energy (E^) such that the extracted value of mf, is 
biased toward larger values because of shape function 
effects [22]. However, this shift cannot be calculated 
model independently. Rather than include a model de- 
pendence, we have multiplied the theory uncertainties 
for the b ^ sj rates by the ratios of the energy differ- 
ence from the endpoint, relative to that for BELLE with 
Ecut = 1-8 GeV.3 

To summarize, we define the combined experimental 
and theoretical error matrix for a given observable to be 



and /o = /i = 1, /2 = 1/4, h = l/(6\/3). Here p are 
the experimental errors, Bi = X^^^' or ^(2") are the co- 
efficients of the last computed terms in the perturbation 
series, and A contains the errors discussed earlier. We 
take A = 0.001 for the data used in the fit, except for the 
CLEO and BABAR photon moments, where we multiply 
A by 1.3 and 1.5, respectively, to account for the increase 
in shape function effects as one limits the allowed region 
of the photon spectrum. 

We stress that there is no unique way to estimate the- 
oretical uncertainties to a given expression. Thus, while 
we believe that our estimates are reasonable, it is cer- 
tainly not the only possible way to estimate the theory 
uncertainties (e.g., taking the theory correlation to be 
identical to the experimental correlations is just an edu- 
cated guess). 



B. Experimental correlations 

Some of the experimental correlation matrices have 
negative eigenvalues. In some cases, these are at the 
level of round-off errors. To avoid these negative val- 
ues, we have added 0.01 to the diagonal entries for the 
correlation matrices for the BABAR and CLEO lepton 
moments, and the DELPHI hadron moments. 

The correlation matrix for the BABAR hadronic mo- 
ments [4] contains negative eigenvalues which are much 
larger than any round-off uncertainties. This persists 
even if only every second value of the cut is used, as 
advocated in [4], so we are forced to add 0.05 to the di- 
agonal entries of the correlation matrix for the BABAR 
hadron moments to make the eigenvalues positive. Note 
that the correlation matrix can have negative eigenvalues 
only if the probability distribution can take on negative 
values. 

The preliminary correlation matrix for the BELLE lep- 
ton and hadron moments was used in the fit [40]. 



no sum OTii,j. 



(16) 



where i and j denote observables, Cij is the experimental 
correlation matrix, and 



for the n'th hadron moment , 

for the n'th lepton moment , 

cTi = \l{<Jr? + {Afn{mB/2Yf + {B,/2f 

for the n'th photon moment , (17) 



C. Constraints on parameters 

Even though there are many more observables than 

there are parameters, the fit does not provide strong con- 
straints on the 1/m^ parameters. Thus it is useful to add 
additional information to ensure that the fit converges to 
physically sensible values of the nonperturbative param- 
eters. Thus, as in Ref. [1] we add to the contribution 



A par am 



_fo, \m<mi, 

\[\{0)\-mlY/M^, \{0)\>ml, 

(18) 

where (m^, M^) are both quantities of order Aqcd, and 
(O) are the matrix elements of any of the 1 / m'^ opera- 
tors in the fit. This way we do not prejudice (O) to have 

3 The B ^ Xs-y photon spectrum also receives contributions of any particular value in the range 1(0)1 < 777,3. ^jjg 
order ml /ml, which are negligible corrections for our analysis. we take = = 500 MeV. We checked in Ref. [1] 
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FIG. 1: Fit results for |Vc6| and mt in the ISexp, ISno, kinEXP, kinNO, and kinuG schemes defined in Eq. (1) and in the 
traditional pole scheme. The red doted, blue dashed, and black solid ellipses denote the results at tree level, order Us, and 
a^Po, respectively, corresponding to Ax^ ~ 1- 



Scheme 


2 

'-^theory 




\Vcb\x 10^ 


mj^ [GeV] 


rrib — rric [GeV] 


mc{mc) [GeV] 


Ai [GeV^] 


ISexp 


yes 


50.9/86 


41.4 ± 0.6 


4.68 ± 0.03 


3.41 ±0.01 


1.07 ±0.04 


-0.27 ±0.04 


kiuEXP 


yes 


52.6/86 


41.2 ± 0.6 ± 0.1 


4.70 ± 0.03 ± 0.03 


3.40 ±0.01 ±0.01 


1.09 ±0.03 ±0.03 


-0.19 ±0.04 ±0.04 


ISexp 


no 


148.4/86 


41.5 ±0.3 


4.69 ± 0.02 


3.39 ±0.01 


1.09 ±0.02 


-0.31 ± 0.03 


kiuEXP 


no 


238.8/86 


41.1 ± 0.3 ± 0.7 


4.74 ±0.01 ±0.11 


3.36 ±0.01 ±0.04 


1.15 ±0.01 ±0.11 


-0.33 ±0.03 ±0.11 



TABLE I: Fit results for \Vch\, nit, nic and Ai in the IS'exp and kiuEXP schemes. Our fits in the kinetic scheme use fi^, but 
the result is converted to Ai to help comparison. The first two lines are the fit results including our estimates of the theoretical 
errors, the lower two lines are setting these to zero. The second error for the kiuEXP scheme is the shift due to changing /i from 
1 to 1.5 GeV. 



that the results for \ Vcb\ and rub are insensitive to varying 
between 500 MeV and 1 GeV. The data are sufficient 
to constrain the I /ml operators in the sense that they 
can be consistently fit with reasonable values, but they 
are not determined with any useful precision. The data 
can be fit without including Xparami but then some of the 
1/ml parameters are not of natural size, with values of 
order 0.5 GeV'^. Including Xparam gives a fit with reason- 
able values of the parameters, of order 0.1 GcV"^. The 
contribution of Xparam is rather small, of order 0.1 — 0.2, 
so that Xparam "^^^s uot drive the fit. This shows that 
there are some very flat directions in parameter space 
which are stabilized by including Xparam- have shown 
our final results for Vet and with and without includ- 
ing Xparam i^^ ^^^^^ final rcsults do not depend 



significantly on whether or not Xparam is included. 

Note that the fit performed by the BABAR Collabo- 
ration included the half-integer hadronic moments. We 
have checked that including these moments still leaves 
some i/ml parameters with values larger than natural 
size. We have chosen to not include these moments in 
the fit since they have large theoretical uncertainties. 

VI. FIT RESULTS AND DISCUSSION 

The fit result for \Vcb\ and nif, in the five mass schemes 
defined in Eq. (1) and in the traditional pole scheme are 
shown in Fig. 1. The fit results are shown at tree level, 
order as, and order a^/?o. The kinetic scheme results 
are obtained using /ij etc. in the fit, and then converting 
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FIG. 2: Measurements (blue squares: BABAR [25, 31], green triangles: CLEO [29, 34], black dots: BELLE [28, 32]) and fit 
results for the hadron invariant mass (left) and the lepton energy moments (right) as functions of the lepton energy cut, Ecut- 
For the hadron moments mj denotes {m\), while for the lepton moments BR is branching ratio. Ml is first moment, and M2 
and M3 are the second and third central moments, respectively. The red (dark) shaded regions show the fit error, while the 
yellow (light) shaded regions are our estimates of the theoretical uncertainties from the A terms in Eq. (17). The A term for 
{mx) and {nix) is three times larger than for (m^) and {m%), because of the worse expansion for the non-integer moments. 



the results back to Ai etc. for easier comparison with the 
other schemes. One can see that the I^exp and kinExp 
schemes have better convergence than the pole scheme. 
The main fit results in the IS'exp and kinExp scheme are 
given in Table I. The quoted \Vcb\ values include elec- 
tromagnetic radiative corrections, that reduce** \Vcb\ by 
VQEB = f + (acm/7r) ln(miy//i) ~ 1.007. The remarkable 
agreement between the fit results shows that the main 
difference in the fits is not which short distance b quark 
mass is used, but whether nib — TTfic is or is not expanded 
in terms of HQET matrix elements. 

The uncertainties for the IS'exp and kiiiExp schemes, 
which eliminate TOc, are smaller than for the ISno and 
kiuNO schemes, which use rfTc{mb). This is contrary to 
the claims made in [9], but is not unexpected, since the 
former schemes have only one parameter at leading order 
in l/nib, while the latter schemes have two such param- 
eters. While not expanding in l/rric gives slightly larger 
errors than expanding, the consistency of the central val- 
ues between the two methods shows that one can use the 
l/rric expansion for inclusive B decays. 

One can clearly see that using the kinetic mass for 



rric (the kinuG scheme) does not reduce the uncertainties 
compared to the IS'exp and kiuExp schemes. Also, as 
is now well known, the pole scheme does not work as 
well in inclusive calculations as the schemes which use a 
short distance mass. Thus, in the remainder of this work 
we will present results in the ISexp and in the kiuExp 
schemes. We have carried out the fits in 6 additional 
schemes, including the PS and MS schemes. All of the 
schemes give reasonable fits, but only the PS scheme with 
nib — ''Tic expanded in HQET gives rise to similarly small 
uncertainties as ISexp and kiuExp- 

The charm quark mass enters into the computation, 
and we can extract the value of nic from our fit. The value 
of nib — nic, which is free of the order Aqcd renormalon 
ambiguity, is (in the ISexp scheme) 

TOfc -mc = 3.41 ± 0.01 GeV. (19) 

We can convert this result to the MS mass of the charm 
quark, 

mdmc) = 0.90 ±0.04 GeV, 

mc(mc) = 1.07 ± 0.04 GeV , (20) 

where the two results depend on whether the perturba- 
tive conversion factor is reexpanded or not.^ The reason 



* In the preprint version of this paper and in Refs. [1, 10] the 
inverse of this factor was used erroneously, which enhanced | Vcb I ■ 
We thank O. BuchmuUer for pointing this out. 



^ I.e., the difference between dividing by 1 + aia^ -|- a2a^ and 
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FIG. 3: Measurements and fit results for tlie liadron invariant mass and the lepton energy moments, setting all theory errors 
to zero in the fit. (See the caption for Fig. 2.) The yellow (light) shaded band gives the estimated theoretical uncertainty, as in 
Fig 2. It is not included in the fit, but it can help to decide the significance of any differences between theory and experiment. 



for the large difference between the two results is that 
perturbative corrections are large at the scale Trie- Tak- 
ing the average of the two rric values, and adding half the 
difference between them as an additional error gives 

mc(mc) = 0.99 ±0.1 GeV. (21) 

The difference between the rric values is a nice illustration 
that one should avoid using perturbation theory at a low 
scale, if at all possible. The kinuc scheme uses pertur- 
bation theory at a scale below rric, and suffers from the 
same problem. 

Next, we compare how well the theory can reproduce 
the experimental measurements, focusing on the cut de- 
pendence of individual moments. The results for the 
hadronic moments and the leptonic moments are shown 
in Fig. 2 in the IS'exp scheme. (The DELPHI and CDF 
results are included in the fits, but are not shown, as they 
correspond to Ecut =0.). The red (dark) shaded band 
is the uncertainty due to the errors on the fit parameter. 
The width of the yellow (light) shaded band is the the- 
oretical uncertainty due to higher order nonperturbative 
effects not included in the computation [the A term in 
Eq. (17)]. Within the uncertainties, the OPE predictions 
for all these moments agree well with the data. As we 
explained before, the moments {m^) and {m\) were not 



multiplying by 1 — aias + {a\ — 02)0^. Only the larger value in 
Eq. (20) has been shown in Table I. 



included in the fit. The yellow bands shown for (mx) and 
{rrix) use A = 0.003 as an estimate of the uncertainty, 
a factor of three larger than for the integer moments, 
because of the worse theoretical expansion discussed in 
Sec. inc. 

The agreement between the theory and experiment for 
the third lepton moment is better than our estimate of 
the theoretical uncertainty. This might be an indication 
that we overestimate the theoretical uncertainty for this 
moment. 

The for the fit shows that the theory provides an 
excellent description of the data. In the IS'exp scheme, 
we get = 50.9 for ly = 86 degrees of freedom, so x^/t^ = 
0.59. The standard deviation for x^/^ is = 0.17, 

so that x^/^ = 0-59 is about two standard deviations 
below the mean value of 1.01. This is some evidence 
that the theoretical errors have been overestimated. To 
study the effect of the theoretical uncertainties, we also 
perform fits with all theoretical uncertainties set to zero. 
This fit gives x^ — 148.4 for the IS'exp scheme, and 
X^ = 238.8 for the kiuExp scheme. The resulting fits 
still agree well with the experimental data, as can be 
seen from Fig. 3. The fit results with no theory error are 
given in the lower half of Table I. The values x^/^ — l-'72 
for the ISexp scheme and x^/^^ = 2.72 for the kinsxp 
scheme are significantly greater than one, which is some 
evidence that there are higher order theoretical effects 
which have not been included. 

The calculations in the kinuG scheme [9] were used by 
the BABAR Collaboration [4], to perform a fit to its own 
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FIG. 4: Fit results for |Vc6| and mj in the kinEXP and kinuG 
schemes using fib ~ 1 GeV (blue and black) and using fib ~ 
1.5 GeV (green and yellow), fic for the kiuuG scheme has 
been kept fixed at 1 GeV. The regions correspond to Ax^ = 1 
(black and yellow) and 4 (blue and green). The upper plots 
includes theory errors in the fit, and the lower plot does not. 

data. While we agree with the results of Ref. [9] for the 
lepton energy moments, we are unable to reproduce their 
results for the hadronic invariant mass moments. One 
should also note that Ref. [9] (i) uses as = 0.22 for the 
lepton moments, and as = 0.3 for the hadron moments 
(ii) includes the a^Po corrections (which are known for 
both the lepton and integer hadron moments) only in the 
lepton moments, but not in the hadron moments. 

The kiuExp and kiuNo schemes depend on a choice for 
fih. In the kinuG scheme there is an additional depen- 
dence on and there is no reason why the theoreti- 
cal predictions should be expanded using m^'"(/if,) and 
toJ^™(/Xc) defined at the same scale {fib = Mc)j since all 
that is required is that each fi should be small, so one has 
to choose both fit and fic- To illustrate the sensitivity to 
the choice of fi,b,c we show a fit in Fig. 4 varying fib from 1 
to 1.5 GeV keeping fic = I GeV fixed. Clearly there is sig- 
nificant dependence in the kinetic schemes with respect 
to changes in /z, and this should be included as an ad- 
ditional uncertainty for that scheme. We have included 
this scale uncertainty in Table I. The kinetic schemes 
use perturbation theory at a low scale fib^c, and so are 
sensitive to precisely how these corrections are included, 
as was the case for Wic{mc)- The PS scheme is much 
less sensitive to the value of fi. In Fig. 5, we show the 
variation with fi in the PS scheme. Note that one advan- 
tage of the IS' scheme is that it does not depend on any 
factorization scale parameter fi. 

Reference [9] quotes smaller theoretical errors than the 
estimates used here, as can be seen from the plots in 



FIG. 5: Fit results for |Vci,| and mb in the PS scheme using 
fi = 2 GeV (blue and black) and using fi — 1.5 GeV (green 
and yellow) The regions correspond to Ax^ = 1 (black and 
yellow) and 4 (blue and green). The upper plots includes 
theory errors in the fit, and the lower plot does not. 

Ref. [4]. We do not believe that this optimistic estimate 
of the theoretical uncertainty is justified. 

Figure 6 shows the results for \ Vcb \ and nib in the ISexp 
and kiuExp schemes with and without including our es- 
timate of the theoretical uncertainties. This plot also 
shows for comparison nib = 4.69 ± 0.03 GeV extracted 
by Hoang [41] from sum rules [42, 43] that fit to the 
BB system near threshold, and the PDG 2004 value [24] 
\Vcb\ = (42.0 ± 1.1 ± 1.8) X 10-3 from exclusive decays. 
Hoang's determination of ml^ is independent of the cur- 
rent determination, and the agreement is remarkable. 
The PDG 2004 value for \Vcb\ from exclusive decays is 
also independent of our determination from inclusive de- 
cays. 

In summary, we find the following fit results: 

\Vcb\ = (41.4±0.6±0.U^) X 10-^ 

= (4.68 ± 0.03) GeV, (22) 

from the ISexp fit including theory errors, where the 
first error is the uncertainty from the fit, and the second 
error (for \Vcb\) is due to the uncertainty in the average 
B lifetime. From the ISexp fit with no theory errors, 
and using the PDG method of scaling the uncertainties 
so that is unity, we obtain 

\Vcb\ = (41.5 ±0.4 ±0.1,^) X 10-^ 

mf = (4.69 ± 0.02) GeV. (23) 

The increase in \Vcb\ compared to Ref. [1] is largely due 
to an increase in the experimental values for the semilep- 
tonic B decay rate since two years ago. 
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0.8 
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B{B^XJi>) i?o(0,0) 


1.0523 ± 0.0001 


1.1216 ±0.0003 


1.246 ± 0.001 


1.466 ±0.001 


1.880 ± 0.003 


B(B ^ XJu,Ecut) Ro{Ecut,0) 



TABLE II: The ratio of the semileptonic branching ratio B{B Xciv) to the semileptonic branching ratio with a lepton energy 
cut Ecut obtained using the 15exp scheme, including theoretical uncertainties. 
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FIG. 6; Fit results for |Vc6| and rah in the I^exp and kinEXP 
schemes. The upper plots include our estimate of the theo- 
retical errors, the lower ones set them to zero. The black and 
blue regions are the fit results (Ax'^ = 1,4), and the yellow 
and green regions are the fit results (Ax^ = 1,4) omitting 
Xparam in Eq. (18). We have also shown a red point given by 
combining Hoang's determination of and the PDG 2004 
value of I Veil from exclusive decays. 



The l^EXP fit (including theoretical uncertainties) also 
gives 



1^^/ ' = (2.49 ±0.02) X 10"" GeV. 



' ^QED 



(24) 



The ratio of the semileptonic branching ratio with no en- 
ergy cut to that with an energy cut is given Table IL 
The semileptonic branching ratio obtained from the fit 
(including theoretical uncertainties) is B{B — > Xc£iy) = 
0.105 ± 0.003. Note that this number depends on the 
PDG 2004 value (corrected for B contamination, 
see Eq. (3)) of 0.105 ± 0.003, and the BABAR branch- 
ing ratio measurements with an energy cut, which give 



a higher value of 0.107 ± 0.002, when converted to the 
branching ratio using Table IL 

Another useful quantity is the C parameter, needed for 
the B A"s7 rate [44], which is defined to be 



C 



T{B Xc 



cb\ 



T{B ^ Xulv) 



0.58 ±0.01. (25) 



The value of C depends on the (unknown) matrix ele- 
ment of four-quark operators, which enter the B Xu(-v 
rate at order 1/m^ (but not B — > XJm). These ab- 
sorb a logarithmic divergence in the 1/m^ corrections 
to the B Xc^v rate in the formal limit uric — > 0. 
(For rric ^ Aqcdj four-quark operators only enter the 
B —t Xclv rate at order ag.) The four-quark operator 
matrix element gives an uncertainty in addition to that in 
Eq. (25). An estimate of the four-quark operators' con- 
tribution is obtained by replacing the formally divergent 
term in the B — > XuiP rate, 8{pi/ml) ln{m'^/ml) [15] by 
8(AQCD/mfc)3 ~ 0.01 [45]. 

The above fits give a robust value for \Vcb\ and m^. 
However, we recommend using the error estimate with 
caution. As we have pointed out, the fit seems to indicate 
that the unknown higher order corrections are smaller 
than our theoretical estimate of 0.1%, so that one can 
use Eq. (23). A theoretical uncertainty less than 0.1% is 
very small for a hadronic quantity at the relatively low 
scale of around 5 GeV. It is interesting that the current 
fit shows that the theoretical uncertainties in inclusive B 
decay shape variables are so small. If this is confirmed 
by further comparisons between theory and experiment, 
the uncertainty in Vcb can be reduced still further. 
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